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ABSTRACT 


Recent observations of the spectra of cosmic ray helium, M, LH 
and VH nuclei in the energy range frcm £ 200 MeV/nuc to > 22 BeV/nuc 
are reported. The differential spectra of all of these groups of 
nuclei are found to have a maximum at 300-400 MeV/nuc at sunspot 
minimum. The average exponents on the integral rigidity spectra in 
the range 5 to 50 BV are 1.54 + 0.03 for He nuclei, 1.50 + 0.04 for 
M nuclei, 1.47 + 0.06 for LH nuclei and 1.40 + 0.08 for VK nuclei. 

The spectra of these groups of nuclei are compared and it is found 
that the average He/M, He/LH and He/VH ratios are 16+1, 70+3 
and 200 + 15 respectively. These values are reasonably constant 
frcm the highest down to the lowest energies measured although some 
evidence for a "dip" is present in the 500-1000 MeV/nuc range for both 
the He/LH and He/VH ratios. Solar modulation effects on these ratios 
are discussed and it is concluded that the ratios measured at earth 
are representative of those existing in interstellar space only if 
energy loss processes in interplanetary space are unimportant. 

The influence of interstellar propagation on the spectra and charge 
ratios at low energies arising frcm ionization energy loss and 
nuclear spallation during matter traversal is examined. It is 
found that propagation models that contain a large number of short 
path lengths significantly modify the expected effects of ionization 
energy loss at lew energies. Specifically it is suggested that the 
presently measured charge ratios are consistent with the passage of 
the average cosmic radiation through enough matter to reproduce the 

3 

abundances of the so-called fragmentation nuclei, Li, Be, B and He . 
Two component models are not required to explain our data. Rather 


we feel that a better representation of the situation results fran 
considering a continuous spatial distribution of sources which, along 
with the actual interstellar propagation conditions, leads to a 
particular distribution of matter path lengths. It is pointed out 
that large differences exist in the approaches used in the literature 
to calculate the effects of matter traversal in interstellar space 
at lew energies. These differences play an important role in the 
interpretation of the experimental results. Significant modifications 
of the charge ratios at lew energies can also be obtained by requiring 
that sane of the matter traversal occur in the cosmic ray sources 
themselves during the cosmic ray acceleration process. This may be 
sufficient to produce charge ratios that are essentially flat at 
low energies even in the presence of interstellar ionization loss. 
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INTRODUCTION 

/ 

/ 

Because of their hiqh charae M(6 Z _< 9) nuclei, LH (10 Z £ 14) 
nuclei and especially VH (Z >_ 20) nuclei should he sensitive indicators 
of the effects of ionization energy losses occurring durina the passacre 
of cosmic ray nuclei throucrh matter. A study of the enertjy spectra 
of these nuclei and a canparison with the spectrum of nuclei of 
lower charge but similar charge to mass ratio, such as helium, is 
capable of providing unique information on the propagation and 
possibly the acceleration of the cosmic radiation. 

Until recently the observed numbers of cosmic ray nuclei with 
>_ 20 have not been sufficient to permit a detailed study of their 
energy spectrum. 

This situation has been remedied by contemporary counter (Webber 
and Ormes, 1967 - Paper I) and emulsion (Freier and Waddincrton , 1968) 
studies of the VH spectrum. The agreement between these separate 
measurements is, in general, satisfactory, although both suffer frcm 
a paucity of measurements of the intensity of VH nuclei above 1.5 
BeV/nuc. In the papers of Freier and Waddincrton where their impressive 
results on the VH spectrum were compared with other workers' measure- 
ments of the helium spectrum, the accuracy of the comparison is 
probably as much limited by the accuracy of the deduced helium spectrum 
as by uncertainties in the VH spectrum. 

Observations of VH nuclei at still lower energies than could be 
reached in these balloon experiments (_< 200 MeV/nuc) have also been 
reported from a satellite by Ccmstock et al., (1966, 1969) and a rocket 
(Reames and Fichtel, 1967) . The rocket measurements give intensities that 
appear to be a reasonable extension of the higher energy balloon results, 
however the satellite observations give differential intensity values 

- - 




by a factor of approximately 2 above 300 MeV/'nuc . As 
a consecruence, the satellite point at - 200 MeV/nuc must be treated 
with reserve. This is unfortunate, since the lowest enerny VH 
nuclei measurements are most sensitive to the possible propagation 
and acceleration effects. 

Although VH nuclei cure the most sensitive indicators of the effects 
of ionization energy loss these effects should nevertheless be 
readily observable for M and LH nuclei as well. Measurements on 
these latter groups of nuclei can be made more accurately and to 
lower energies than for VH nuclei. Therefore they provide a useful 
confirmation of the VH data as well as enabling the systematics of 
any variations with charge to be examined. 

We have recently extended our M, LH and VH nuclei measurements 
reported in Paper I, and have obtained 5 new integral intensity 
values above 1 BeV/nuc and extending out to > 22 BeV/nuc. In addition 
the low energy parts of the M, LH and VH spectra have been measured 
more accurately than in Paper I. In this paper we shall compare 
these new measurements on M, LH and VH nuclei with our recent measure- 
ments of the helium spectrum covering the range frcm m ioo MeV/nuc 
to > 24 BeV/nuc (Ormes and Webber, 1968) . Much of this helium data 
has not been published previously. In most instances the helium, M, 

LH and VH nuclei spectra have been measured simultaneously thus 
eliminating the need to correct for the sometimes considerable 
effects of solar modulation. 

INSTRUMENTATION AND BALLOON FLIGHTS 

The data to be discussed here was obtained on balloon flights 
using similar, but ouite distinct, instruments. The details of 


the balloon fliqhts are criven in Table 1. The operation of these 
instruments is as follows: 

Instrument 1. This is basically a Cerenkov-scintillation 

counter telescope . It contains two thin scintillators which form. 

the telescope coincidence, followed by a 15 an dia. UVT. lucite 

Cerenkov detector, subsequently followed by another (ranne) scintillator. 

The first three telescope elements are all pulse height analyzed 

usirn separate 2048 channel analyzers. This instrument, its response 

to charged particles and the method of unfolding the differential 

spectrun from the measured pulse heiqht distribution have been 

described in detail in a previous publication, (Ormes and Webber, 

1968) . It is used to obtain the differential enemy spectra of 

helium, M, LH and VH nuclei from ^ 100 fieV/nuc to 2 BeV/nuc. The 

Cerenkov threshold of 300 MeV/nuc for lucite, and the enemy oorrespond- 
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ing to a range of 8.8 a/am provide calibration points for the low 
energy nuclei. This 4 -element system provided a substantial improve- 
ment in background surpression and charge resolution over the earlier 
two element system upon which the results of Paper 1 were based. 

The lamer area and loncrer flight times provided more than 2 times 
the total area- time factor than in the entire earlier series of 
measurements. 

Instrument II: This is a lame area double scintillation counter 

telescope (diameter of individual elements = 61 a ) . This system 
is used to obtain the intearal intensity of nuclei at gecmaonetic 
latitudes < 40° where only relativistic nuclei of each charne are 
present. The scintillators are optically compensated to obtain 
uniformity in light collection to within + 2%. The technique for 
compensation , and the details of the charge resolution etc. are 





it is described in Paper I 


being reported in a separate publication (von Rosenvinge and Webber, 

1968) . The analysis of flight data using this instrument is 
particularly simple. Because of the good charge resolution and 
low background, each charge stands out as an individual clump of 
particles in the 2 -dimensional matrix of scintillator outputs. 

Instrument ill . This instrument is identical to II except that 
a 1.7 meter long gas Cerenkov detector is inserted between the 
scintillation counter telescope elements. Thus, in addition to 
the integral intensity of particles above the rigidity threshold 
determined by the geomagnetic field at the location of the flight, 
the integral intensity of events above the gas Cerenkov threshold 
is obtained. In the two measurements with this instrument reported 
here, Freon gas was used at two atmospheres and one atmosphere pressure, 
giving effective thresholds of 15 and 22 Bev/nuclecn respectively: 

(for helium these thresholds are 16 and 24 BeV/nuclecn) . A paper 
describing the operation of this system and a more complete discussion 
of the results obtained, is new in preparation. For these studies 
it is relevant to note that, since the output of the gas Cerenkov 
detector is pulse height analyzed, the gas threshold energy can be 
rather clearly defined, particularly for heavy nuclei. 

In the last analysis the ability to study the relative spectra of 
individual charges depends on the degree of resolution of the different 
charges and on the statistical accuracy of the results. This, in turn, 
depends on the amount of "background" present in the pulse height 
matrix and the intrinsic resolution of the detector for adjacent charges. 
A pulse height histogram of all events above the Texas cut-off for the 
charge range 2*10-18 for the two Texas flights with instrument III is 
shown in Figure 1. Similar charge resolution is obtained for the other 
flights including those at Ft. Churchill using instrument I. 
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ATMOSPHERIC CORRECTIONS 

The observed intensities at balloon altitudes have been extrapolated 

2 

to the top of the atmosphere usim absorption mean free paths of 50 a/cm , 

33 g/an 2 , 30 a/an 2 and 16.5 a/cm 2 respectively for Heliun, M, IH and VH 

nuclei as directly measured in air, (Paper I) . It is assumed that these 

mean free paths are enemy independent following the work of Cleahom 

(1967) . The correction for ionization energy loss in the atmosphere 

has been applied separately. It should be noted that the mean depths 

2 

of overlying matter were between 2.5 and 6 g/an for the various flights 
so that atmospheric op r recticn& are not believed to be unduly serious. 

FSSXTS 

"he heliun and VH nuclei results are sumarized in Table II. The 
total nunbers of VH nuclei actually observed are also listed so that 
one may evaluate the statistical errors. 

The differential enemy spectra for these two nuclei are shown 
in Figure 2. The dashed curves represent the anoo th ed spectra (measured 
or deduced) used in the analysis of Freier and Waddimton (1968) . The 
published results of Ccmstock et al. , (1966, 1969) , and Keanes' and Pichtel 
(1967) on VH nuclei at low enemies are also shown. In addi t ion, there 
are measurements on VH nuclei by Bhatia et al. , (1966) and Dumaprasad 
and Reanes (1968) which provide valuable indications of the behaviour 
of the VH spectra. These measurements are in good agr eemen t with the 
more statistically accurate balloon measurements and are not shown in 
Figure 2. There are no major differences in the new spectra we pr e sen t 
here and those used by Freier and Haddington. However, there are 
differences in detail as seen by a comparison of the dotted and solid 
lines in Figure 2. These lead to an appreciable difference in a 
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TABLE II 


OUSNSHY OF HELIUM AND VH NUCLEI 
(Epoch 1966) + 


Etergy 

or Rigidity He Nuclei VH Nuclei 

2 

(particles/in -ster-sec) 


190 - 

240 MeV/nuc 

12.5 + l a 


- 

240 - 

280 MeV/nuc 

10.5 + l a 

(9) 

0.050 + 0.018 

280 - 

440 MeV/nuc 

44.5 + 3 a 

(52) 

0.195 + 0.03 

440 - 

560 MeV/nuc 

29 + 1.5 a 

(30) 

0.146 + 0.03 

560 - 

700 MeV/nuc 

26 + 1.5 a 

(33) 

0.155 + 0.03 

700 - 

1000 MeV/nuc 

40 + 2.5 a 

(51) 

0.230 + 0.036 

> 1000 

MeV/nuc 

119 + 3 a 

(290) 

0.64 + 0.05 

> 15 

BeV/nuc 

4.63 + 0.3 b 

(8) 

.026 + 0.010 

> 22 

BeV/nuc 

2.46 + 0.2 C 

(8) 

0.015 + 0.006 

> 3.20 BV 

126 + 5 

(39) 

0.60 +0.11 

> 4.35 BV 

91.5 + 2.0 

(115) 

0.445 + 0.04 

> 5.00 BV 

79.7 + 1.6 

(163) 

0.39 1 + 0.03 

>11.1 

BV 

25.1 + 0.8 

(67) 

0.13 5 + 0.015 


a. Derived from spectra presented by Qrmes and Webber (1968) . 

b. Effective threshold * 16 BeV/nuc. 

c. Effective threshold = 24 BeV/nuc. 

4#ft. Washington neutron monitor * 2350. 


comparison of the two spectra, particularly at low enemies. 

The differential spectra of M and LH nuclei obtained in this study 
are shown in Fioure 3. The satellite results of Ccmstock et al., (1966, 
1969) , Balasubrahmanyan et al. , (1966) and the rocket observations 
of Reames and Fichtel (1967) are also shown at lav enemies. 

The integral ricridity spectra for helium and VH nuclei are 
shown in Fioure 4, and for M and LH nuclei in Figure 5. Above 5 BV 
rigidity the exponents on these spectra are: helium nuclei = 1.54 + 
0.03, M nuclei = 1.51 + 0.05, IH nuclei = 1.47 + 0.06, VH nuclei = 

1.40 + 0.08 . In deriving these spectra a mass to chame (A/Z) ratio 
of 1.95 is taken for helium nuclei, 2.00 for M and LH nuclei and 
2.15 for VH nuclei. 

The ratios of the differential energy spectra of helium and VH 
nuclei are shown in Figure 6. The earlier results obtained on this 
ratio by Freier and Waddington (1968) are also shown in this figure. 

The peculiar "kink" at 600-800 MeV/nuc. which appears in their 
results and complicates the interpretation of the enemy dependence 
of this ratio is still present to seme extent in our data. Our more 
detailed measurements at higher energies permit a clearer picture of 
the overall variation of this ratio to be obtained, however. 

In Figures 7 and 8 we show the He/M and He/Til ratios obtained in 
the series of measurements. Note that in the case of the He/IH ratio 
there is a "kink" at v 400 MeV/nuc similar to that appearing in 
the He/VH ratio at v 700 MeV/nuc. This "kink" is not apparent in 
the He/M ratio, however. 

SOLAR MODULATION EFFECTS 


The differential spectra and ratios just presen ted represent 



the situation at the top of the atmosphere . Before we can make 
astrophysical deductions frcm the measurements we must attempt to 
correct for the effects of modulation in the solar environment. 

We do not have a completely adequate theory of solar modulation 
that will permit us to demodulate the respective spectra to obtain 
the spectra existina in interstellar space. It is generally 
accepted that the spectra of nuclei with the sane A/Z ratio will 
be modified identically by the modulation process. The nuclei we 
are concerned with here have A/Z ratios which are different in seme 
case and therefore solar modulation effects may distort the 
relative spectra. By simultaneously measurino the spectra of helium 
and heavier nuclei we have eliminated a source of distortion inherent 
in comparisons of these nuclei at different levels of solar modula- 
tion. However, residual solar modulation effects operative at 
the time of measurement cannot be removed by this procedure and 
indeed have not been considered in earlier discussions of the charrre 
ratios. 

Two aspects of the solar modulation are of importance in the 
"demodulation” of the spectra measured at earth to those existing 
in interstellar space. They are: (1) The functional (enemy or rigidity) 
dependence of the modulation, and (2) the magnitude of the residual 
modulation existing at sunspot minimum. The current state of our 
experimental understanding of these aspects has been reviewed most 
recently by Webber (1968) and we shall follow the conclusions reached 
therein. 

The modulation M of a particular nuclear species is taken to be 
the In of the ratio of intensities of that species at times t^ and t 2 
at a given kinetic energy/nucleon . Let us call R ftB (t) the ratio of 




the differential intensities of species A and B at time t also at 
the same kinetic energy/nucleon. The modulation of a given species 
is described in the usual diffusion convection picture by a form 
M = where K(t) is a quantity describing the level of the 

solar modulation (sometimes called the modulation parameter) , D is the 
diffusion coefficient and P is the particle rigidity. The functional 
dependence of the modulation is contained in the diffusion coefficient. 

It is easy to show that velocity dependent solar modulation (eg. 

D % 3 ) will not change the ratio R^g. Measurements indicate that 
above ^ 200 MeV/nuc. the functional dependence of the diffusion coefficient 
is more closely ^ 6P. The rigidity dependent factor will result in 
the ratio R^ changing with the level of solar modulation. We have 
calculated this change for a residual modulation parameter, = 0.6 BV. , 
which is believed to be representative of the value necessary to demodulate 
spectra measured in 1965-66 at earth to those existing in interstellar 
space, (Gloeckler and Jokipii, 1967) . The increase in the He/VH ratio 
that results is 3% at 600 MeV/nuc. , 8% at 200 MeV/nuc. and 15% at 
100 MeV/nuc. The effects are less for the He/M and He/LH ratios 
because of their more nearly identical mass to charge ratios. These 
modifications are comparable to or less than the experimental errors 
in the ratios themselves and we have therefore not adjusted the ratios 
near the earth shown in Figures 6, 7 and 8 for this effect. 

It has been pointed out by Webber (1968) that certain experimental 
data are not in agreement with the predictions of the diffusion- 
convection model for the modulation of cosmic rays. Indeed a recent 
paper by Gleeson and Axford (1968) has investigated this problem and 
develops a modulation picture in which the effects of energy loss apparently 
determine the solar modulation. The predictions of such a model with 
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recrard to the variation of the respective chame ratios are substantially 
different than in the r'iffusion-convection picture w^ere it is assumed that 
we are seeina narticles at earth at the sare enerov at vhich they exist 

in interstellar space. These differences will be more nronounced at lower 

* 

enemies, fccordincr to Gleeson and Oxford (1968) the enemy loss 
corresponding to the residual modulation at sunsnot nininun may be ^ l n 0 
Pe'Vnuc. Particles of 200 be w /nuc observed at earth miaht then corres- 
pond to 300 M eV/nuc particles in interstellar space (ie. the ratios of 
Figures 6, 7 and. 8 would be shifted ^ 100 M eV/nuc to the rinht to be 
representative of interstellar space) . Pronaeational effects on the 
He/T*, He/JJI and. He/VH ratios at lm enemies such as ionization enemy 
loss that could be detectable in interstellar space miaht then be 
unobservable at earth since me would be unable to see the appropriate 
enemy particles. 

At this staae of our understandina of the solar modulation process 
it is not feasible to mabe more detailed calculations based on the 
above hypotheses. It should be obvious that uncertainties in both 
the narmitude and functional dependence of the residual solar modulation 
are considerable and it is conceivable that the modifications of 
the chame ratios at low enemies may not be nealiaihle. mhis is 
certainly true in the case of t he actual spectra themselves where 
residual modulation effects may cmpletely alter the shape of the 
low enemy portion of the spectra. For this reason we shall emphasize 
a study of the He/ M , He/LH and He/VH ratios rather than a discussion 
of the features of the individual spectra themselves. 

For the purpose of further discussion we shall consider the ratios 
presented in Figures 6, 7 and 8 as representative of those existing in 
interstellar space. T,1 e shall new consider briefly their interpretation. 
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The He/M. He/LH and He/VH Ratios and the Propagation of Cosmic Rays 
through Matter in Interstellar Space 

It has been generally argued that heavy cosmic ray nuclei must 

2 

have travelled through an average of 4 g/cm of matter (taken to 

be interstellar hydrogen) (eg. Shapiro and Silberberg, 1968; Beck 

and Yiou, 1968). This argument is based on measurements which show a 

substantial intensity of light nuclei (Li, Be, B) in cosmic rays. These 

light nuclei can be regarded as products of the fragmentation of heavier 

nuclei as they pass through matter. Calculations based on the observed 

2 3 

intensities of other "fragmentation" nuclei such as H and He lead to 

the same conclusion; that cosmic rays in general must have, at some 

2 

time in their lifetime, passed through ^ 4 g/cm of matter. 

The effects of the passage of low energy heavy nuclei through 

2 

the equivalent of 4 g/cm of matter without a concurrent energy 
gain process are profound indeed. The enormous rate of energy 
loss by ionization of these heavy nuclei relative to helium will 
severely depopulate the low energy part of the spectrum and cause 
the initial He/M, He/LH and He/VH ratios to increase markedly at 
low energies. The problem then becomes how to reconcile the charge 
ratios at low energies which do not exhibit a sharp increase and there- 
fore seem to suggest that only a small amount of matter is traversed and 
the measurements of the absolute abundance of fragmentation nuclei which 
imply passage through considerable matter. 

The expected increase in the charge ratios depends importantly on 
the characteristics of the propagation of cosmic rays in interstellar 
space. It also depends on the initial or injection spectra which 
might not be identical. To make the problem tractable it is usually 
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assumed that all charge species Initially have the same spectra. Further- 
more this spectrum Is usually assumed to be a power law In total energy. 

In this way the effects of propagation may be examined seperately. 

The cosmic ray propagation problem has recently been treated 
in a comprehensive manner by Fichtel and Reames (1968) . These authors 
have examined the effects of certain simplified propagation models on 
the various charge ratios. They solve the fundamental transport 
(diffusion) equation of cosmic rays in the galaxy including the effects 
of energy loss by ionization and fragmentation due to nuclear inter- 
actions. They consider certain limiting path length (material) 
distributions based on specific boundary conditions. Three path length 
distributions they consider seem to be appropriate for a comparison 
with the data. All of these distributions will as a first requirement 
reproduce the observed abundance of the "fragmentation" nuclei 
Li, Be and B. 

The first distribution is the simple slab length approximation in 

which all particles of a given energy traverse the same amount of 

matter, a 6 - function approximation. The expected charge 

2 

ratios when the matter thickness ■ 4 g/cm , constant with energy are 

shown as curves in Figures 6, 7 and 8. Here the full effects of 

matter traversal at low energies are observed, the resulting ionization 

loss causing the charge ratios as shown in these figures to increase 

rapidly at low energies. The data is clearly inconsistent with this 

type of path length distribution-indeed only a small fraction of the 

2 

matter traversal (< 1 g/cm ) could actually occur in interstellar space 
on this picture. 

The second type of path length distribution arises from a point 
source diffusing into an infinite medium. This produces a roughly 
gaussian path length distribution for particles of a fixed energy. 


a 
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This smearing of paths lengths tends to destroy the features of the 

charge ratios which depend on the specific energy dependence of the 

nuclear spallation cross sections (cf the He/VH ratio in Figure 6). 

The expected charge ratios when the average of this gaussian 

distribution is 4 g/cm , constant with energy, are shown as curves 

A2 in Figures 6, 7 and 8. The increase in the ratios at low energies 

is not as pronounced because of the existence of an appreciable 

fraction of short path lengths. However, the data is still, generally 

speaking, not in good agreement with such a path length distribution. 

The third path length distribution arises from a uniform spherical 

distribution of sources with the sun at its center. Cosmic rays then 

diffuse throughout this region. The expected charge ratios when 
2 

R 0 /X(R 0 = radius of region, X - mean free path for diffusion) = 9.3 
2 

g/cm , constant with energy are shown as curves in Figures 6, 7 
and 8. Here the distribution of path lengths is strongly weighted towards 
short path lengths, hence the increase in the ratios at low energies 
is least of all. In fact this model gives a reasonable fit to the 
measurements for all three charge ratios, He/M, He/LH and He/VH. 

At this point it is certainly relevant to comment on the experi- 
mental errors inherent in the measurements of the charge ratios at low 
energies. In view of the large differences obtained by various 
observers (eg. values of the He/M ratio range from 17 to > 25,* values 


*For example Comstock et al. (1969) from data taken aboard the 
OGO-I satellite find a He/M ratio actually decreasing markedly at 
low energies, whereas Fan et al. (1968) on the Imp-Ill satellite present 
an He/M ratio almost constant with energy. Moreover neither of these 
sets of ratios are consistent with those obtained on a rocket by Reames 
and Fichtel (1967) or reported by Hagge et al. (196E) (see Fig. 7 for 
data points) . 
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of the He/LH ratio from 60-100, at energies ** 200 KeV/nuc) it is clear 
that these uncertainties are substantial. This definitely complicates 
the interpretation of the results. 

The third model discussed above is probably the most realistic 
from the point of view of our current understanding of the cosmic ray 
origin-propagation problem and the magnetic field structure of the 
galaxy (eg. Davis, 1962; Ginsburg and Syrovatskii, 1964) although 
it does net represent a completely acceptable picture of this 
situation by any means. This model serves to illustrate the importance 
of considering a distribution of sources, however, and clearly represents 
a model which provides a path length distribution skewed to short path 
lengths. As such it minimizes the increase in charge ratids at low 
energies due to ionization energy losses. 

This path length distribution will predict a decrease in the 
expected L/M ratio at low energies because of the preponderence of 
short path lengths. The curves presented by Fichtel and Reames (1968) 
show that the L/M ratio should decrease from ^ 0.29 at 4 BeV/nuc to 
0.18 at 100 MeV/nuc. We believe that the data on the L/M ratio 
presented by von Rosenvinge et al. (1969) is not inconsistant with 
such a decrease. Thus we believe that all of the current data on 
the spectra and relative abundance of helium and heavier nuclei can 
be explained in terms cf Fichtel and Reames' model 3 or some 
other model of source distribution and propagation characteristics 
that results in a similar path length distribution. 

The path length distribution arising from model 3 of Fichtel 


and Reames is in fact quite similar to the exponential path length 
distribution model suggested by Cowsik et al. (1967). The Cowslk model 
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ls equivalent to considering propagation in a bounded medium and will 
contain an even larger fraction of short path lengths. Obviously the 
charge ratios at low energies will be even less effected by ioniza- 
tion energy loss; however, the variation of the L/M ratio to be 
expected at low energies is now too extreme to be consistent with 
the data on this ratio presented by von Rosenvinge et al. (1969) . 

This model must therefore be rejected on this basis. 

Comstock (1969) has examined the variation of the charge ratios 
at low energies using OGO-’i data only and Waddington and Freier (1968) 
have examined the variation of the He/VH ratio over a wide range of 
energies. All of these workers consider that a sharp increase in 
these ratios should be observed at low energies as a result of ioniza- 
tion energy loss. The data presented by these workers does not 
exhibit the expected increase (in broad agreement with our results). 

In both instances the authors have suggeeted, however, that the reason 
for this is the presence of a more or less well defined "second" 
component of cosmic rays at low energies. It is postulated that this 
"second" component traverses very little matter-hence no upturn in the 
charge ratios at low energies - and also has the proper chemical abundances 
and energy spectra to fit smoothly onto the known spectra for the high 
energy component. 

He would like to draw attention to several differences between 
the calculations in the above papers and those of Fichtel and Reames 
regarding the propagation of cosmic rays in Interstellar space and 
to consider alternatives to the approaches which lead to the two 
component models. The first point deals with the procedure used 
by both Comstock (1969) and Haddington and Freier (196$ for the 
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calculatior of the effects of Ionization energy loss. Both of these 
papers consider a 5-function, or slab approximation, for the path 
length distribution. One would expect, then, that this calculation 
would be similar to model 1 developed by Fichtel and Reames (1968). 

There are important differences in the two sets of calculations, 
hdwever. Fichtel and Reames solve a 3 dimensional transport 
equation, in the case of model 1 with the somewhat artificial condition 
of a 6-function path length distribution. The calculation of the 
change in spectrum due to ionization energy loss made by Comstock 
(1969) and Waddington and Freier (1968) is a one dimensional slab 
length approximation. Waddington and Freier carry their calculation 
no further than this while Comstock goes on to integrate the solutions 
over one dimensional Gaussian and exponential path length distributions. 

The difference between the one dimensional and three dimensional 
slab approximations is quickly realized by considering an example. 

Energy loss by ionziation is a dissipative process. A simpler energy 
loss process to consider, for the sake of an example, is that in a 
(conservative) electric field. Suppose therefore that there is 
a one-dimensional electric field (corresponding to a positive 
potential V) between the source and the observer. Following the 
reasoning of the slab approximation, we find that J OBSERVER 
(E-ZeV/A) - JgouRcg( E )- 11118 18 Liouvllle '» theorem for one- 
d intension. On the other hand, Liouville's theorem for 3-d imens ions 
tells us that 

J OBSERVER (E-ZeV/A) J SOURCE <E) 

— - - ' - constant 

P (E-ZeV/A) p vE) 
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in this case the observer spectrum is multiplied by an additional 

2 2 

factor [p (E-ZeV/p (E)] <1. The effects of energy loss on the spectra 

2 

are amplified by this factor. At low energies, since p ^ E, the spectral 
change calculated using the one dimensional and 3 dimensional approaches 
will differ by a factor ^ E. A comparison of curves A^ and M in 
Figure 6, which differ by a factor ^ 1/E at low energy, shows that 
this additional factor is also apparently contained in the dissipative 
calculations as well. 

| Comstock then goes on to integrate the one-dimensional solutions 

over a distribution of path lengths. There are then two problems to 
this approach: 1) Given a particular geometrical model, what is the 

proper distribution of path lengths? 2) Is the 3-dlmensional Liouville's 
theorem (dissipative case) in fact satisfied? Comstock in no way 
addresses himself to either problem but merely considers two hypothetical 
path length distributions. The advantage of the Flchtel and Reames 
method is that Liouville's theorem in three dimensions is automatically 
satisfied (it is the starting point) and the machinery for obtaining 
the correct distribution of path lengths is also set up directly. 

In order to justify a one dimensional calculation one has to 
believe that during the propagation of cosmic rays they are rigidly 
tied to the galactic magnetic field limes with essentially no diffusion 
Davis (1962) has studied this possibility and in fact presents 
arguments in favor of a highly "ordered" motion of this type. It is 
therefore possible that the one dimensional approach used by Comstock 
(1969) and Haddington and Freler (1968) is more nearly correct. 

The limitation in the approach of these authors is therefore not so 
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much in the type of propagation considered but that they do not 
consider realistic distributions of path lengths arising from a 
specific source distribution-magnetic field configuration. They have 
considered very specialized situations. The results of Fichtel 
and Reames (1968) show that the changes in charge ratios at low 
energies depend importantly on the characteristics of the propaga- 
tion and of the source distribution itself. 

He believe therefore that the difficulties Comstock and 
Haddington and Freier try to explain by a two component model are 
more readily understood in terms of a continuous path length 
distribution. It should be emphasized that a continuous spatial 
distribution of sources may be as important in determining this 
path length distribution as the actual propagation conditions 
themselves . 

Before completing our discussion we should like to consider a 
further alternative to understanding the behavior of the various 
charge ratios at low energies. 

If the ratios presented in Figures 6-8 are indeed found to 
remain essentially flat or even turn down at the very lowest energies 
(< 100 MeV/nuc) then apparently this behavior cannot be completely 
explained by path length distribution considerations either in a one 
dimensional or a three dimensional approach. In this case two 
distinct component source models such as suggested by Comstock (1969) 
and Haddington and Freier (1968) are certainly an alternative 
explanation, however a simplier alternative is to retain a basic 
one component source model and allow for some traversal of matter 
in the source region Itself. Under certain conditions, placing 
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a fraction of the matter in the cosmic ray source is equivalent to 
providing for a preferential acceleration of heavier nuclei at lower 
energies. This arises when (dE/dt)/ accelerat ^ is only slightly 
larger than (dE/dt)/ ionization i oss * or the highest charges accelerated. 
In this situation VH nuclei are not accelerated as rapidly to 
higher energies as are He nuclei. Energy loss by ionization does 
not depopulate the low energy end of the VH spectrum either. The net 
effect is an excess o£ VH nuclei relative to He at energies just above 
the injection energy as compared to the situation at high energies. 

The mechanics of such a situation for a Fermi accelerated total 
energy spectrum in the source have been discussed by Ramadurai (1967), 
who has illustrated the modification brought about. One of us (T. T. 
von Rosenvinge) is presently considering models in which some but not 
all of the material traversal occurs in the sources during acceleration 
and effectively obliterates the effects of ionization loss. 

SUMMARY AND CONCLUSIONS 

The results of this paper may be summarized as follows: 

(1) We find the He/M, He/LH and He/VH ratios are essentially 
constant with energy in the range > 100 MeV/nuc up to > 22 BeV/nuc. A 
"dip" that may be significant occurs in the He/LH and He/V.I ratios at 
500-1000 MeV/nuc. The average ratios above 100 MeV/nuc are 16+1, 

70 + 3 and 200 + 15 respectively. Large uncertainties (discrepancies), 
previously not pointed out, exist in the data on the charge ratios at 
lower energies as obtained on satellites. 

(2) It is argued that solar modulation effects will not 
appreciably alter these charge ratios above a few hundred MeV/nuc. 
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At lower energies the Interstellar ratios could be considerably altered 
If the solar modulation can be described In terms of an energy loss. 

This might make an upturn that is present In interstellar space 
unobservable at earth. 

(3) The influence of Interstellar propagation on the changes 

in the relative spectra at low energies arising from lonzlation energy 
loss during matter traversal has been examined. Propagation models that 
contain a large number of relatively short path lengths significantly 

modify the expected effects of ionization energy loss at low energies. 

i 

| Specifically it is found that the presently measured constancy of 

charge ratios is consistent with the passage of the average radiation 

I 

I through considerable matter. At the present time "two" source 

f 

| models are not required to explain this data and it seems to us 

§ 

§ more appropriate to think in terms of a spatial distribution of 

sources. 

(4) It is pointed out that large differences exist in the 
approaches used in the literature to calculate the effects of 
matter traversal at low energies. These differences play an 
important role in our understanding of the propagation of cosmic 
rays and in our Interpretation of the experimental results. 

(5) Significant modifications of the charge ratios at low 
energies can be obtained by requiring that some of the matter 
traversal occur in the cosmic ray sources during the cosmic ray 
acceleration process. This may be sufficient to produce charge 
ratios what are essentially flat at low energies even in the 
presence of Interstellar ionization loss. 
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FIGURE CAPTIONS 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Charge distribution of nuclei from Z ■ 10 to 28 observed in 
Texas in May- June, 1967. 

Differential energy spectra of helium and VH nuclei (Epoch, 

1966) . Our data points shown as open and solid diamonds 
respectively. VH spectrum multiplied by 100. Dashed curves 
are the spectra of these nuclei used in the analysis of Freier 
and Waddington 1968. N represents the average of 1964-1965 
VH measurements of Reames and Fichtel (1967) . C the 
measurement of Comstock et al. (1969). 

Differential energy spectra of M and LH nuclei (Epoch, 1966). 

Our data points shown as open and solid diamonds respectively 

Points labeled C are Chicago 0G0-I results (Comstock et al. 
o 

1969); Cj are Chicago IMP-III results (Fan et al. 1968); 

N„ are NASA rocket results (Reames and Fichtel, 1967); N are 
NASA satellite results (Balasubrahmanyan et al. 1966); and 
N q are NASA 0G0-I results (Hagge et al., 1968). 

The integral rigidity spectra of helium and VH nuclei (Epoch, 

1966). VH spectrum is multiplied by 100. Helium and VH data 
from this experiment shown as solid rectangles, VH data from 
Paper I as solid squares with error bars. Instrimen tally deter- 
mined energies below 1 BeV/nuc adjusted to rigidities using a 
mass/charge ratio m 2.0 for He and 2.15 for VH nuclei. Freier 
and Waddington (1968)VH measurements shown as open circles. 

The integral rigidity spectra of M and LH nuclei (Epoch, 1966) . 

M nuclei data shown as solid circles, LH nuclei data as open circles. 
Instrumentally determined energies below 1 BeV/nuc adjusted to 
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rigidities using a mass/charge ■ 2.0 for both M and LH nuclei. 

Figure 6 The ratio of helium to VH nuclei as a function of energy. 

Our data points shown as solid diamonds. The explanation 

of the data symbols is as follows: C^, Chicago OGO-I results 

(Comstock et al. 1969), N^, NASA rocket results (Reames and 

Fichtel, 1967) smoothed curve through Minnesota results 

of Freier and Waddington (1968). Theoretical dashed curves 

Aj_, A 2 , and A^ are from the work of Fichtel and Reames (1968) 

as described in the text. Dashed curve M is the calculation 

2 

for a slab length * 4 g/cm (equal total energy spectra) by 
Waddington and Freier (1968). . 

Figure 7 The ratio of helium to M nuclei as a function of energy. Our 
<~ata points shown as solid diamonds. The explanation of the 
lata symbols is as follows; Cj, Chicago IMP III results (Fan 
et al. 1968); C Q Chicago OGO-I results (Comstock et al. 1969), 
N^, NASA rocket results (Reames and Fichtel, 1967), N, NASA 
satellite results (Balasubrahmanyan et al. 1966) and N Q , NASA 
OGO-I results (Hagge et al. 1968). Theoretical dashed curves 
Ap A 2 and A 3 are from the work of Fichtel and Reames (1968) 
are described in the text. 

Figure 8 The ratio of helium to LH nuclei as a function of energy. 

Our data points shown as solid diamonds. The explanation of 
data symbols and theoretical ratios is the same as in Figure 7. 
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